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Figure 1 (b) 
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Figure 2 (a) 
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Figure 2 (b) 
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Abstract 

The Sachs- Wolfe temperature fluctuations produced by primordial density perturba- 
tions are proportional to the potential field (j), which is a weighted integral over the density 
field 6. Because of the central limit theorem, (j) can be approximately Gaussian even 
when 6 is non-Gaussian. Using the Wold representation for non-Gaussian density fields, 
£(r) = //(|r — r'|)A(r')d 3 r', we find conditions on A and / for which (j) must have a 
Gaussian one-point distribution, while 8 can be non-Gaussian. Sufficient (but not neces- 
sary) conditions are that the density field have a power spectrum (which determines /) 
of P(k) oc k n , with —2 < n < +1, and that A(r) be non-Gaussian with no long-range 
correlations. Thus, there is an infinite set of non-Gaussian density fields which produce a 
nearly Gaussian one-point distribution for the Sachs- Wolfe effect. 



The COBE observations of Sachs- Wolfe fluctuations in the cosmic microwave back- 
ground place severe constraints on possible models for large-scale structure (Smoot et al. 
1992; Wright et al. 1992). One of the unresolved questions regarding the primordial density 
fluctuations which gave rise to large-scale structure is whether these density fluctuations 
were Gaussian, i.e., whether the N-point distribution of primordial density fluctuations is 
a multivariate Gaussian distribution. Gaussian models have the virtue of simplicity, and 
a great deal is known about their properties (Bardeen, Bond, Kaiser, &: Szalay 1986); on 
the other hand, a large number of non-Gaussian models have also been proposed. Conse- 
quently, many recent papers have addressed the issue of looking for a non-Gaussian signal 
in the COBE data (Kung 1993; Srednicki 1993; Luo &: Schramm 1993; Luo 1994a,b; Hin- 
shaw et al. 1994; Smoot et al. 1994), as well as predicting the microwave fluctuations to 
be expected in non-Gaussian models (Bennett &: Rhie 1993; Coulson, et al. 1994; Pen &: 
Spergel 1994; Moessner, Perivolaropoulos, &: Brandenberger 1994). 

However, Sachs- Wolfe observations cannot provide a decisive test of whether or not 
the initial density field is Gaussian; in this paper, we derive conditions under which a non- 
Gaussian density field can produce a Gaussian one-point temperature distribution. We 
mention first the well-known possibility that the density field could be non-Gaussian on 
small scales but Gaussian on the length scales probed by COBE. Such density fields are pro- 
duced, for example, by a Poisson distribution of objects; the density field is non-Gaussian 
when smoothed on scales smaller than the mean separation of objects, but Gaussian when 
smoothed on scales larger than this mean separation (Scherrer &: Bertschinger 1991). 

We can, however, make a stronger assertion: the potential field can have a Gaussian 
one-point distribution when smoothed on a given length scale even when the density field 
is non-Gaussian on this same scale, even in the limit where the smoothing length goes to 0. 
The Sachs- Wolfe fluctuations probed by COBE are directly proportional to the potential 
field 0(r) on the surface of last scattering. This potential field is related to the field of 
density fluctuations 6(r) via 

V 2 0(r) = ATrGa 2 p6(r), (1) 

where a is the cosmological scale factor and p is the mean density. [Although the Sachs- 
Wolfe effect is inherently relativistic, our use of the Newtonian expressions for the potential 



will give the correct answer for the problem of interest here. This is justified in the 
Appendix]. In integral form, we can write this relationship as: 

m = -G*j^. (2) 

Thus, the potential field is a weighted integral over the density field. The central limit 
theorem indicates that this integral can be driven to a Gaussian distribution even when 
the distribution of 8 is non-Gaussian. A similar phenomenon has been noted in the case 
of the velocity field, which can also be expressed as an integral over the density field; 
the distribution of velocities can be Gaussian even when the underlying density field is 
non-Gaussian (Scherrer 1992). 

The first step in our argument is to express the density field in terms of the Wold 
representation (Scargle 1981; Peebles 1983). Any density field can be represented as the 
convolution of a spherically symmetric function / with a white-noise function A: 

*(r) = |/(|r-r'|)A(rVV. (3) 

The function A is a white noise function in the sense that the two-point correlation function 
is zero: (A(i"i)A(r2)} = for i"i ^ ^2- However, A can (and in general will) have non- 
zero higher-order correlations. The proof that any density field can be represented in the 
form of equation (3) is trivial. The Fourier decomposition of 6(r) gives us a set of Fourier 
components £(k). If we write each £(k) as £(k) = R{\s.)e ^ \ where R and 9 are real, 
then the Fourier transform of R gives / (which must be spherically symmetric because R 
is real) and the Fourier transform of e gives A (which has a white noise power spectrum, 
since |e ( )| 2 = 1, but which may have non-zero higher-order correlations). In a sense, 
the Wold decomposition is the most natural representation for non-Gaussian models. The 
phase correlations in the model are contained entirely in the function A, while the power 
spectrum is given entirely by the function /. For example, the phase correlation statistic of 
Scherrer, Melott &: Shandarin (1991) involves calculating the properties of the field A(r). 
A number of authors have used the Wold representation to generate non-Gaussian fields 
(Peebles 1983; Lucchin & Matarrese 1988; Moscardini et al. 1991). 



We will assume that the density field on the scales probed by COBE has a power 
spectrum P(k) oc k n . The Fourier components of the density and potential field are given 

by 

6(k) oc k n/2 A(k), (4) 



anc 



^(Jfe) ex k n / 2 ~ 2 A(k), (5) 

where A(k) is the Fourier transform of A(r) which appears in equation (3). To derive 
an expression for /(|r — r'|) in the Wold representation of the density field, we Fourier 
transform k n ' 2 to obtain 

£(r)oc f A(r')\r-r'\- n / 2 - 3 d 3 r'. (6) 

[This Fourier transform of a power-law in k is discussed by Peebles (1980). Strictly speak- 
ing, convergence of the Fourier transform here requires n/2 > —3 and smoothing at small 
scales; equation (6) is then the correct expression for |r — r'| greater than the smoothing 
scale. Also, equation (6) fails for the white noise case n = 0; in that case we simply have 
6(r) = A(r).] Equation (1) or (2) tells us that (j) can be written as 

0(r)cx f Air'^r-r']- 71 / 2 - 1 ^^. (7) 

[Again, this expression is valid only for n/2 > —1, which is satisfied by the COBE results 
(Wright, et al. 1994; Gorski et al. 1994). The case n = +4 corresponds to a white-noise 
potential distribution, for which (j)(r) oc A(r).] 

We will assume that the density field and potential field are smoothed with some 
window function of radius R, where R can be arbitrarily small. Note that since smoothing is 
a convolution, we can simply smooth A(r) with this window function to give the correctly- 
smoothed density and potential fields. Now we have to make some assumption regarding 
the nature of A. For now, we will assume that A has a non-Gaussian one-point distribution 
p(A), when smoothed over the length scale R, but that all higher-order correlations for 
r > R vanish. In a sense, this is a minimally non-Gaussian model; such models (in which A 
was used as the density field) were investigated numerically by Messina et al. (1990). This 



also describes a special class of the seed models investigated by Scherrer &: Bertschinger 

(1991) and Scherrer (1992), namely, seed models in which all of the seeds have identical 

mass and are uncorrelated. In this case £>(A) is just a discrete probability distribution and 

/ is the accretion pattern around each seed. 

We can now apply the central limit theorem for the sum of weighted independent 

random variables (Feller 1971, p. 264). Consider a random variable X with a known 

distribution, having zero mean and unit variance. Now construct a set of independent 

random variables Xk given by Xk = &kX . If the random variable Zjy is the sum of the 

first N Xfc's: 

N 

Zn = 2_^Xk, (8) 



and sn is the variance of Zj^: 



k=i 



N 



s 



N = }] a li 
k=l 



(9) 



then Zn/sn has a Gaussian distribution in the limit N — > oo if both 

s 2 N — > oo, (10) 

and 

a N /s N -> 0. (11) 

By smoothing on a length scale R, we have effectively discretized the potential and 
density fields. Think of the A field as divided into discrete regions of volume ~ R 3 . Then 
equations (6) and (7) can be written as sums, where we translate to the origin: 

oo 

S(0)oc^A(r fc )|r fc r"/ 2 - 3 i? 3 , (12) 

jfe=i 

and 

oo 

^(Ojoc^AMIr*!-"/ 2 - 1 ^ 3 . (13) 

jfe=i 
We choose an order of summation such that we first sum over all regions at a given radius 

|r|, then move on to the next larger radius. This gives |rj.| oc k 1 ' 3 , so a for the j term in 

each sum is 

^■oc^aJ-"/ 8 - 1 , (14) 



for the density field, and 

^occtaJ-/ 6 - 1 / 3 , (15) 

for the potential field, where a\ is the variance of the A field. Then the quantity s 2 N is 

just 

N 

^occriVV"/ 8 - 2 , (16) 



for the density field, and 



k=i 



N 
2 



<r 2 A J2 k ~ n/3 ~ 2/3 ( 17 ) 

k = l 

for the potential field. Consider first the case of a Zel'dovich power spectrum [n = +1). 
For the density field, we note that s 2 N oc 5^jfe=i k~ 7 ' 3 , which converges to a finite value 
as N — > oo. Thus, the central limit theorem will not apply to the density field in this 
case. For the potential field, we have s 2 N oc Y2] e=1 k~ 1 , so that sjy — > oo as N — > oo. 
Furthermore, <tn oc TV -1 ' 2 so <tn/sn —> as N — > oo. Thus, the central limit theorem is 
satisfied, and the potential field will have a Gaussian distribution. 

We can ask the more general question: what are the limits on n for which (j) must be 
Gaussian (independent of the one-point distribution of A), while 6 can be non-Gaussian? 
We note that the sum ^ fc k s converges for s < —1 and diverges otherwise. Then equation 
(17) tells us that for the potential field, sjy — > oo for N — > oo as long asn< +1, while for 
n < +1, <tn/sn ~ 1/v-ZV for large N, so equation (11) is also satisfied. Conversely, if we 
wish to violate equation (10) for the density field, then equation (16) tells us that n > —3. 
These results make intuitive sense. For —3 < n < +1, the integral over the power spectrum 
which gives (£ 2 ) is dominated by large k (small length scales), while the corresponding 
integral for ((j) 2 ) is dominated by small k (large length scales), so it is reasonable that the 
central limit theorem works for (j) but not for 6 in these cases. Since equation (7) is valid 
only for n > —2, we can conclude that for a power spectrum of the form P(k) oc k n , with 
— 2 < n < +1, there is an infinite class of models for which the density field can be highly 
non-Gaussian, while the potential field has a nearly Gaussian one-point distribution. This 
range includes the Zel'dovich power spectrum and is consistent with part of the allowed 
range given by COBE (Wright et al. 1994; Gorski et al. 1994). Note that o"a must be finite 



for our argument to be valid, although no physically-realistic models have been proposed 
which produce an infinite variance for the density and potential fields. 

To illustrate this result, we have generated two sets of density and potential fields with 
different power spectra, but the same non-Gaussian distribution for A(r). We begin by 
setting down the field A(r) on a 128 3 lattice with a gamma function one-point distribution 
p(A) oc Ae , and no correlations between the values of A on different lattice sites. We 
Fourier transform A(r) to obtain A(k), and then we take £(k) = fc 1 ' 2 A(k) and (j>(k) = 
fc" 3 / 2 A(k) for Fig. 1, and 6(k) = k 1 ^) and 0(k) = jfe-^k) for Fig. 2. (Multiplicative 
constants are irrelevant for our purposes). We then inverse Fourier transform 6(k) and </>(k) 
to obtain the density field 6(r) and the corresponding potential field (j)(r). To eliminate 
lattice effects, we smooth both of these fields with a Gaussian window function W(r) = 
e~ r2 l 2R \ with R = 3 cell lengths. 

Our figures give the one-point distribution of the density and potential fields for 
P(k) oc k n with n = 1 (Fig. 1) and n = 2 (Fig. 2). The distributions of both density fields 
(Figs, la and 2a) closely resemble the original gamma function. This is not suprising, since 
in both cases the form for /(|r — r'|) in the Wold representation (equation 3) is sharply 
peaked near zero. The corresponding distributions for the potential fields are shown in Figs, 
lb and 2b. In the Zel'dovich case (Fig. lb) the distribution of the potential field retains 
a slight imprint of the original gamma distribution, but it is nearly indistinguishable from 
a Gaussian, confirming our argument that the potential field will be driven to a Gaussian 
in this case. Fig. 2b does not satisfy our conditions for the central limit theorem to apply, 
and the distribution of the potential in this case is obviously non-Gaussian. 

Although we have derived our results for the case where A is uncorrelated, our ar- 
gument generalizes to the case where A has short-range correlations, but is uncorrelated 
over distances greater than some finite length (Ibragimov &: Linnik 1971). Suppose that A 
has non-zero correlations for separations less than some length ro , but is uncorrelated over 
distances greater than ro- Then we can group terms in equation (13) into regions with 
some size r ~^> ro, and let tpj be the contribution to the sum in equation (13) from the j 
region: (j) = J^ ■ tpj. There will be correlations between the tpj's for adjacent regions, but, 
since r ^> ro, there will be no correlations between regions which are not adjacent. Now 



we break the sum in equation (13) into two pieces, (j) = (j)\ + (j>2, dividing our regions in 
checkerboard fashion, so that no two regions in the sum for (j)\ are adjacent, and no two 
regions in the sum for $2 are adjacent. The central limit theorem applies to both (j)\ and 
(j>2] both must have a Gaussian distribution. Now we fix the size of the regions for (j>2, but 
increase the size of the regions for (j)\ to arbitrarily large size; the result is that (j) % <j>i, so 
(j) has a Gaussian distribution. 

Note that this argument does not require smoothing: it applies even to the case where 
the smoothing length is much smaller than the maximum length over which A is correlated. 
[In the opposite case, the argument is trivial.] This is a special case of a more general result 
given by Ibragimov &: Linnik (1971), who show that the central limit theorem applies if 
the correlations go to zero sufficiently rapidly with increasing separation. Thus, it is likely 
that our results can be generalized further, to sufficiently weakly correlated A, although 
what "sufficiently weakly" means in this context is unclear. 

Our results have several limitations. As we have noted, they do not apply to cases 
where A has long-range correlations or infinite variance. Of course, our results cannot 
be completely general, since it is obviously possible to choose the potential field to be 
non-Gaussian by construction. Furthermore, our results apply only to models in which the 
dominant cosmic microwave fluctuations are proportional to the potential at the surface of 
last scattering, not to models in which secondary effects are important. Finally, we have 
shown only that the one-point distribution function for the potential field is Gaussian, 
although it may be possible to generalize our results to the full N-point distribution func- 
tion. [Non-Gaussian models with a Gaussian one-point distribution are certainly possible, 
but somewhat contrived (Feller 1971; Scherrer, Melott, &: Shandarin 1991).] Within these 
limitations, however, our results do indicate that there is an infinite set of models for which 
the density field has a non-Gaussian distribution, but the potential field has a Gaussian 
one-point distribution, even when smoothed on arbitrarily small scales. Thus, even if the 
cosmic microwave background turns out to have a Gaussian one-point distribution on a 
given scale, this observation by itself cannot prove that the underlying density field is 
Gaussian on that scale. 

R.J.S. was supported by NASA (NAGW-2589) and by the Department of Energy 
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APPENDIX 

In calculating the temperature fluctuations due to the Sachs- Wolfe effect, we must 
evaluate the gravitational potential on scales larger than the horizon size at decoupling. 
Since many quantities, most notably the density contrast (6p/p), depend on the choice of 
gauge under these conditions (Bardeen 1980; Kodama &: Sasaki 1984; Mukhanov, Feldman, 
&: Brandenberger 1992), we must be careful to insure that our results are gauge-invariant. 
We address the gauge issue here. 

The Newtonian gravitational potential has a natural gauge-invariant generalization, 
usually denoted by an upper case $ [called $h by Bardeen (1980)] to distinguish it from the 
usual Newtonian potential (0). The gauge-invariant Sachs- Wolfe formula linearly relates 
temperature fluctuations to $, so our results hold true if $ has the same statistics as (j). We 
will show this to be true. However, a possible source of confusion arises. The gauge invari- 
ant generalization of the density perturbation is not unique. In fact the gauge-invariant 
density fluctuation variable (called e g by Bardeen 1980) which is used by Mukhanov et al. 
(1992), reduces (in flat space) to —2$ on scales larger than the horizon during the matter 
dominated era, which seems to imply that there is no difference between the density and 
potential field, hence apparently nullifying our conclusions. 

One way to see that our conclusions are correct is to use a different generalization 
of the density fluctuation called e m . This is the density fluctuation relative to the space- 
like hypersurface which represents the local matter rest frame everywhere. This is the 
most natural variable for looking at matter perturbations. With this variable the usual 
Newtonian-like relation 

V 2 $ = AnGpe m (Al) 

holds for all scales, including suprahorizon scales. Thus no problem arises as the relation 
between e m and $ does not change with scale. The apparent contradiction stems from 
the fact that when one constructs the variable e g , one cancels out the gauge dependence 



by adding something which contains the potential itself. On subhorizon scales all gauge- 
invariant definitions of the density contrast converge to the same answer. On suprahorizon 
scales, however, the relationship between the two definitions mentioned here is 



e 9 



K 

1 



Qa 2 H 2 



ld$ / K 

H~dt + V~ a 2 E" ' 



(A2) 



where K is the scalar curvature, a is the scale factor, is the fraction of the critical 
density, H is the Hubble parameter, and t is the time since the big bang. In flat space, 
(K = 0), the potential $ for perfect fluids is constant, regardless of whether the universe 
is dominated by matter or radiation. Hence eq. (A2) reduces to 

e g - e m ~ -2$. (A3) 

Thus the reason the variable e g is ~ —2$ on suprahorizon scales is that the —2$ has been 
explicitly added on during the process of trying to cancel out the gauge dependence of 
6p/p. 

With this distinction between different definitions clear, we now argue that the choice 
of variable has nothing to do with the argument put forth in the main part of the paper. On 
subhorizon scales, all gauge-invariant generalizations will reduce to the usual Newtonian 
6p/ p. Thus, if the density fluctuations responsible for the present subhorizon structures 
are indeed non-Gaussian in character, the present potential field $ must have the near 
Gaussian distribution under the conditions specified in this paper. If we extrapolate $ 
back to the time of recombination, $ will still exhibit this Gaussian behavior, since the 
potential is approximately constant in time. It is then irrelevant which density fluctuation 
variable we use to describe the density perturbations at recombination. The potential 
$, and hence the temperature fluctuations, will have the Gaussian distribution under the 
conditions outlined in the text. 
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FIGURE CAPTIONS 



Fig. 1: (a) The distribution of densities in the model given by the Wold representation: 6(r) = 
/ /(|r — r'|)A(r')rf 3 r', where A(r) is an uncorrelated field with a gamma distribution, 
and / is chosen to give a Zel'dovich power spectrum [P(k) oc k] for the density field. 
Solid curve is a Gaussian distribution with the same mean and variance, (b) The 
distribution of potentials for the same model. Solid curve is a Gaussian distribution 
with the same mean and variance. 

Fig. 2: As Fig. 1, for the power spectrum P(k) oc k 2 : (a) Distribution of densities, (b) 
Distribution of potentials. 
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